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Abstract 

In this work the Feynman-Schwinger representation for the two-body Greens 
function is studied. After having given a brief introduction to the formalism, we 
report on the first calculations based on this formalism. In order to demonstrate 
the validity of the method, we consider the static case where the mass of one of the 
particles becomes very large. We show that the heavy particle follows a classical 
trajectory and we find a good agreement with the Klein-Gordon result. 



1 Theory 

Recently, the Feynman-Schwinger representation (FSR) was presented |T[] as a new co- 
variant formalism to calculate the relativistic two-body Greens function. It was shown 
that the FSR sums up all ladder and crossed diagrams and that it has the correct static 
limit if one of the masses of the two particles becomes very large. Furthermore it was 
argued that the formulation was well suited for essentially nonperturbative gauge the- 
ories such as QCD, since the formalism can be set up in a gauge invariant way and 
the possibility exists to include a vacuum condensate in the interaction kernel via the 
cumulant expansion. While the formalism neglects all valence particle loops, one has the 
possibility to include all self-energy and vertex-correction graphs as well. In this paper 
however, we will concentrate on the static limit. 

In order to demonstrate the formalism we consider the 3 -theory for two charged 
particles xi an d X2 with masses mi and vri2 and charges g\ and 5% interacting through 
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exchange of a third, neutral particle with mass fi. The Euclidean action for this theory 
is: 



s 



E 



d 4 x 



= d 4 x 



xl (m 2 -dl + g i( p) Xi + \<t> (V ~ dj) 



(1) 



where summation over the index % is implied. The Greens function is defined as the 
transition probability from the initial state = x\{ x i)X2{x 2 ) to the final state \I/f = 
Xi(yi)Xa(l/2): 

G^!,^^!,^) ex | V X iV X 2V(t) V t % e" 5 * 

(detA 1 (y 1 ,x 1 )A 2 (y 2 ,x 2 ))-^ A^i, xJA^fa, x 2 ) e"^ d4 ^ 2 -^ (2) 



Next, in order to get the sum of all generalized ladder graphs, we neglect the determinant 
in @. This corresponds to neglecting all the xi~ an d X2-loops and is often called the 
'quenched approximation'. We now wish to rewrite (|2|) in such a way that we can 
perform the integration over the remaining field <fi as well. To this end we exploit the 
Feynman-Schwinger representation for A -1 : 



A \y,x) 
ds lim 



x 



N y 



ds (y \ e - s{m2 - 9 * +9 ^ 

2N „JV-1 

lb _L 

ds J (Vz) xy e~ m2s -^ Jo ^ {r)dT - gs Jo t( z ( T )) dT 

Integrating over the 0-field is now straightforward and yields: 
G(y 1 ,y 2 \x 1 ,x 2 ) oc 



zn=v 



ZQ=X 



(3) 



— mjsi— rrij 



^ ds 2 J (V Zl ) Xiyi (Vz 2 ) X2y2 e 

(4) 

In case we leave out the self-energy contribution and the vertex-corrections, the Wilson 
loop (W,},) is simply given by: 



exp 



exp 



9i92SiS 2 



N 2 



N 

a f i (zij + — ^ 2 ,fc — ^ 2 ,fc~i 



j,fe=i 

-1 /■! 



9i92SiS 2 f f A (^i(ri) - 2 2 (r 2 )) dridr 2 



where A is the free scalar two-point function: 



A(z) 



47T 2 U 



■K^\z\) 



(5) 



(6) 



This formulation for the Greens function has the great advantage that it is essentially 
a quantummechanical one in which the two valence particles interact via a nonlocal 
interaction. One can show that (|4]) and (|5]) sum up all ladder and crossed diagrams. 
Furthermore, if m 2 — > oo then G obeys the following Klein-Gordon equation: 



(^-m 2 -ff^(|r|)) G(r) = *«(r) 



(7) 



with 



V(\r\) 



dt A 



+ t 2 



which reduces to the instantaneous Coulomb interaction for the case n —>■ 0. 
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2 Results 

For a Euclidean theory it is known that for large times T = T^MA^Mh^^yMzEMl 
the contributions to G are dominated by the lowest lying states and that they fall off 
exponentially, G = J2iCiexp(—EiT). Hence, we can obtain the ground state energy by 
observing that: 

Eq = - lim -^logG = - lim ( -^-G) IG (8) 

The time derivative with respect to T can be done explicitely. 

We project the Greens function on a complete set of total and angular momentum 
states |P, l,m), but due to the nonlocal interaction (|5|) this does not imply that the 
degrees of freedom that are associated with the generators of the symmetries (P, L 2 , 
L z ), can be integrated out. 

For all our calculations we used the Metropolis Monte Carlo algorithm to perform 
the integrations over s\, S2 and all the coordinates. The value of N in (|3|) and @ that 
we needed to get reasonably stable results was typically 25, so that effectively we had 
to perform a 200-dimensional integral. Convergence was usually reached after 3 ■ 10 8 
points which took approximately 30 hours of CPU-time on our fastest workstation. In 
order to demonstrate that it is indeed possible to determine the actual groundstate 
within this formalism, we have investigated the case of one light and one heavy particle 
interacting through exchange of a massless third particle (mi = 1, mi = 100 and fi = 0). 
For the strength of the coupling we took g 2 = —g\g2 = 4500. The distributions of the 
coordinates of particle 2 are confined to a very narrow region around their classical values, 
demonstrating that the heavy particle indeed follows a classical trajectory. Particle 1 
shows much more quantum behavior; its distributions have finite widths. The normalized 
distributions of S\ and S2 are shown in figure [l|. Note the difference in scales between 
both distributions. 

In figure |2| we show the average of 5 runs as a function of the number of Monte Carlo 
points. The exact result for m2 — > oo, obtained by solving the Klein-Gordon equation 
(0), is Eq = 100.55 and is indicated by the arrow on the right hand side of the figure. The 
FSR-result fluctuates around this value and the average over the last 4 • 10 8 Monte Carlo 
points is Eq = 100.57 ± 0.02. This clearly demonstrates the feasibility of the method. 

We have also studied in detail the situation where mi = m<i = 1 and /i = 0. This case 
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O -theory 

m,= l, m 2 =100, T=500, N=25, g 2 =4.5*10 3 




MC-steps 

Figure 2: Bound-state mass in the static limit as a function of the number of Monte 
Carlo points. T is given in units of the inverse total mass. 

is well suited to test the ladder approximation since the solution of the Bethe-Salpeter 
equation in the ladder approximation is known exactly 0, [|, §]. For values of the 
coupling constant such that there is a binding energy of roughly 10% of the mass of the 
constituents, we find that the full results lie substantially below the ladder predictions. 

In conclusion, we have demonstrated that the FSR is well suited for going beyond the 
ladder theory in the study of composite systems in field theory. In particular, ground 
state properties such as the binding energies and bound-state wave functions can readily 
be determined with this method. 
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